Transport in Graphene superimposed by a moving Electrical Superlattice Potential 



(N 

O 

(N 

(N 



Jiirgen Dictel^ and Hagen Klcincrt^'^ 

'institut fiir Theoretische Physik, Freie Universitat Berlin, Arnimallee I4, D-14195 Berlin, Germany 
^ICRANeT, Piazzale della Repubblica 1, 10 -65122, Pescara, Italy 
(Dated: Received July 3, 2012) 

We calculate dc-conductivities of ballistic graphene undulated by a overlying moving unidirectional 
electrical superlattice (SL) potential whose SL-velocity is smaller than the electron velocity. We 
obtain no dependence of the conductivity on the velocity along the direction of the superlattice 
wavevector. In the orthogonal direction however, the dependence is strong on the velocity especially 
at voltages where a new Dirac point emerges for zero velocity. It is shown that the infinite graphene 
system can serve as an ideal motion detector at potentials where the first new Dirac point emerges. 
There the conductivity is zero at vanishing SL velocities and jumps to infinity when the SL starts 
moving. For finite systems at voltages where the number of new Dirac points is of the order of the 
ratio of the electron velocity by the SL-velocity, the modifications to the conductivity of a moving 
SL is at least of similar magnitude as the conductivity of the stagnant SL. 

PACS numbers: 72.80.Vp, 73.21. Cd, 73.22.Pr 
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I. INTRODUCTION 

The electrical conductivity in suspended graphene 
samples show high mobilities where ballistic transport 
is seen for samples up to the micron length . Due to 
the quasi-relativistic behavior of its electrons, graphene 
has a density of states proportional to the electronic en- 
ergy which is zero at the neutrality point. As a con- 
sequence, this leads in ballistic graphene to the phe- 
nomenon that the conductivity shows a universal finite 
behavior Q whose precise value is still under debate @. 
It seems now that the universal conductivity in a wide 
range graphene sample with highly doped leads has the 
value a = Ae'^/irh (y, |7[, where in a system with vanish- 
ing small doped leads, it is ct = e^7r/2/i Q. Numerically, 
these two values are quite close to each other. A small 
perturbation of the chemical potential of the graphene 
sample may be caused by applying an external gate volt- 
age, this conductivity can change drastically, due the now 
finite density of states at the Fermi-energy. For an infi- 
nite large ballistic system, it becomes even infinite. Such 
an extreme sensitivity of the neutral graphene system 
on the environment parameters makes it attractive as a 
building block for nano-detectors. It was experimentally 
shown that graphene is a good chemical sensor which is 
able to detect the dc-response changes due to the adsorp- 
tion of even single gas molecules on its surface @. This 
high sensitivity is mainly due to the intrinsic low-noise 
properties of graphene. A more general review of possible 
graphene sensors can be found in Ref. [lol 

Here we consider a ballistic graphene sample with an 
overlying slowly moving unidirectional electrical SL. We 
calculate the longitudinal conductivities along and or- 
thogonal to the SL wavevector as a response of a small 
external dc-field. This system is considered as a possi- 
ble model for a graphcnc-bascd nanomechanical motion 
detector. 

In the direction orthogonal to the wavevector of the 
SL we obtain, especially at SL voltages where new Dirac 



points emerges in the non-moving SL, a high sensitivity of 
the conductivity values on the SL-motion. In the parallel 
direction our approximation produces no dependence on 
the SL-velocity. 

Graphene under the effect of a moving SL can be real- 
ized for example by placing periodically patterned gate 
electrodes on either a moving underlying substrate or on 
a rested substrate where now the individual gate elec- 
trodes are activated appropriately with time such that an 
effective moving SL is simulated. More directly, the ex- 
perimental realization could be also carried out by using 
the coupling of the graphene sheet to the electrical field 
of a surface acoustic wave on a piezoelectric substrate 
[ll| or to a charged moving membrane with ripples. 

It was recently shown explicitly for graphene that new 
Dirac points in the energy spectrum can be opened by 
imposing a non-moving SL on the graphene lattice |12| - 
[T^ . This leads to unusual conductivity properties in such 
systems [l5l - [l8| . These new Dirac points are accompa- 
nied with new energy valleys. Due to the technical com- 
plications in handling transport in a moving SL we will 
consider at first the transport contributions of the inner- 
valleys near the K and K' points in Sect. Ill, then those 
of the outer-valleys in Sect. IV. Note, that such a sep- 
aration is not useful for the non-moving SL as will be 
seen in Sect. IV. We start in Sect. II by reconsidering 
first the lowest-band eigenvalues and eigenfunctions for 
the non- moving SL. 

Wc discuss here the most simple representation of a 
SL being a symmetric two-step Kronig-Penney poten- 
tial with a superlattice potential V{x) — Vx{x) where 
x{x) = sg[sin(27rx/(i)] (cf. Fig. 1). The function sg[a;] 
is the sign of x, and d is the wavelength of the SL. In 
the continuum approximation, the graphene Hamiltonian 
under consideration near the Dirac point K is given by 
Hy^ — hvp{axdx/i + crydy/i) + V{x + Vst) [1]. Here Wi? is 
the Fermi velocity and ax,y are the Pauli matrices, while 
Vs is the velocity of the moving SL. 
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II. LOWEST BAND EIGENVALUES AND 
EIGENFUNCTIONS 

In the following, we solve the eigenvalue equation 
H,j^u^''{r') = eu''=(r') for a non-moving SL {vg = 0) 
by using the transfer matrix method [H, [l^ . For the en- 
ergy dispersion in the lowest band we restrict ourselves 
to the lowest-lying oscillatory regime \es\d/hvp ^ V,aQ 
and obtain 1201 



(1) 



with CteA^) = {[(^ 

sin[ao]e*"°/ao: <5o = 



- vix))/hvF]^ - fc2}i/2d/2, r = 

ao/V where V = Vd/2hvp. The 
Bloch momentum in a;-direction is restricted to —ir/d < 
kx < Tr/d. The parameter s = 1 denotes the conduction 
band and s = — 1 the valence band. 

The corresponding lowest-band cigenfunctions are 
(.T, y) = e*'^!'^u^= (x), in the fundamental zone < a; < 
d and for Vg = reduces to u°(a;) = A(a;)u^(0) with 
A(x) = Xo{x)@{d/2 - x) + Xd/2{x)Xo{d/2)e{x - d/2), 
where 



Xxo {x) = cos 



'afXx)2{x - .To)' 



sm 



a^^ {x)2{x—xq) 



d 



E+ 1 iM. 

(2) 

Here E is the unit matrix and M = kyU^ + i[es ~ 
V{x)]a2/hvF- u°(0) is given in the oscillatory region 
\es\d/hvF < V, ao by 



u°(0) 



cos(ao) sin(ao) 




.■sin (gp) 



Vk„d 



(3) 



where iVu in ([3|) denotes a normalization factor. From 
([1]) we obtain an oscillatory behavior of the lowest band 
eigenvalues as a function of ky (cf. Fig. 1). New Dirac 
points emerge at k = for 1/ G Ntt. These are shifted 
along the y-axis in k-space for increasing V. Note 
that the lowest band energy values beyond the oscilla- 
tory regime with momenta ky ^ {V/HvpY scale like 
|e,| ~ hvF\ky\ Bill. 

In the following, we discuss the transport contribu- 
tions of electrons in the inner-energy valleys where ky <C 
V/hvp and the outer- valleys where fio ^ 1 separately. 



III. INNER- VALLEY TRANSPORT 
CONTRIBUTIONS 

In the inner- valley regime ky <^ V/hvp, the lowest- 
band cigenfunctions Uj(a;) for the non-moving system 
stated above are given by 



u°(x) 



1 



kyV* 



hvpkyT' 



(4) 



where Nu in ([4]) denotes a normalization fac- 
tor. The phase factor 4'±{^) is given by 
4>±{x,t) — exp[i5'J° (x, with S\{x,t) = 
Tihj^ dx'sgn[V{x')]aeSx')/{d/2) - iest for = 
chosen such that solves simultaneously the corre- 
sponding time-dependent Schrodinger equation (TSE). 
From (|4|) wc deduce the remarkable observation that the 
inner-valley electrons do not backscatter at the potential 
steps. This phenomenon is well known for ordinary 
Dirac-fermions as Klein-paradox. 

In the following, we use the inner- valley approximation 
ao{x) « V{l-kl/2) with ky = hvpq/V in (g]) that is 
a good approximation of the overall oscillatory behavior 
of the energy dispersion in Fig. 1. Similar approximations 
will also be used when solving the TSE for Vg below. 
Finally we note that the missing of the kx , ky dependence 
in the vector part {^1, 1)"^ of both spinor components in 
^ is due to the inner- valley restriction ky <^ [V/hvpY- 

We obtain from ([T|) that an entire set of 2[t//7r] -I- 1 
Dirac points exists near K where [x\ is the lowest integer 
number smaller than x. By using the inner-valley ap- 
proximation discussed above, these new Dirac points are 
located at /c^d « 2[2V{V - im)^''^ with n = l... [V/i:] 
and fc° = (restricting ourselves to positive ky). The 
linearized energy spectrum around these Dirac points is 
given by e" =shvF[k1 +r^j(fcj^ — kyY]^/"^ where the effec- 
tive y-velocity coefhcient is given by r„ « 2[1 — mi/V] 
for n = 1, . . . , [F/tt], and Fq = sin(V")/V' for the cen- 
tral valley. The magnitude olq for ky ~ ky is given 
by Tm/V. Below, we shall also need the fcy-momentum 
spacings between the right and left-energy crest and the 
Dirac point. The spacing for the right crest is given 
by Ak'j^-^d « 7r[y/2(l/ - 7r?i)]i/2 for n = l,...,[V/7r] 
and Aky'^ — Aky'^ for the left crest positions where 

n = 1 . . . [V/tt] — 1. For the central crest distance we ob- 
tain AfcO'^d = Afc^'^d ^ /r^U^ A] , ,T^„T>2 

and Afcf /-l'^ = fef/^l - AfcO-^ 
the inner-valley formula with ky ^ 1 considered in this 
subsection is valid for the valleys 1 — im/V <^ 1 with 
n 7^ and also the central valley n = 0. 

In the following we solve the TSE ihdtUg" {x , t) = 
Hy^u^' {x,t) with the initial condition u^=(a;, 0) = u°(j:) 
for t = in the oscillatory regime by using the above 
approximations. Note that by using the characteristic- 
method we can solve the TSE without approximation for 
ky = 0. This leads again to (|4]) where now is Vg- 
dependent. Instead of doing this explicitly, we can gen- 
eralize this procedure to any non-zero ky <^ (V/hvp)'^ by 
the Hamilton- Jacobi Ansatz 



y ^-^..y ...Ak\;'-^dr/[iki;'-^dr+AroV'] 

y — "-1/ ^ ^"'y ■ Fin^^lly "we note that 



-hvpsg[V{x + Vst)]^{dxS'^Y + kl + V{x + Vst) 

(5) 

with the boundary condition that S'J'(x,0) — ^^-(a;, 0). 
Due to the local uniformity of V{x + Vst) in position and 
time we obtain local uniform solutions of That this 
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kyd/2 

FIG. 1: (Color online) Lowest Bloch band energy spectrum 
for fci = as a function of dimensionless momenta ky for 
various SL potential strengths V (the full lowest band energy 
spectrum can obtained by using its mirror symmetry with 
respect to the x,y-axis). Here we used the transfer matrix 
method [l^ . [T9| . Inset shows a graphene layer with an over- 
lying moving SL in x-direction. 



approach leads to a TSE solution in the oscillatory regime 
is due to the fact that the general solution can be written 
as 

IfcxI jg--st[ft"FSVfe|+^+V(x+^),t)]gifc,a; 

(6) 



(x, t) ^y^cLs 



in the inner- valley regime fc^. The complex vari- 

ables as_fc^ are local uniform functions in the (x, t)-plane. 
We will show below that as,fc^ is non-zero for only two 
special fca;-values which moreover fulfil the inner-valley 
regime condition k^<^ k^. 

We now solve ([5| by using a generalized characteris- 
tic method for the Hamilton- Jacobi equation that is well 
known in the semi-classical approach to quantum me- 
chanics [2l| . This is based on the one-particle mechanical 
trajectory of a relativistic particle and anti-particle in a 
step-potential. The calculation is outlined in App. A. 

After some manipulation we obtain the result 



S^' with 
and 



Sl-^=±[A±a^*+v:t)+B±ax*Tv*pt)+C±txix*Tv*Ft)], 

(7) 

where ^{x) — dx'x{x'), x* ~ x — [v* — Vs)t and 



V 



1^ 



Vs ± Vf 

V 

{vs ± Vf) 



Vs 
VF 



Vs ± Vf 

k^Vs 



,C± 



±Vk'yZ±, 



zvf Vs ± Vf 



(8) 



--VF 



4)' 



We restrict here the solution of (O to small velocities 



< 



Next we calculate the dc-response in the moving SL 
system. This is done in the gauge A = — cE(t — to)6(t — 
to) assuming to < in general. Since dii{t) does not 
depend on to for t 3> we set immediately to = 0. 
The total Hamiltonian in the contiimum approximation 
is then given by Ha = -ffug + hvF{e/ c){axAx + UyAy). 
The corresponding TSE-solution which we expand to 
first order in A and assume it to satisfy the initial 
condition u^(t = 0) = u° is denoted by u^. From 
this solution we obtain the conductivity in the i-th di- 
rection by da = liriiEi^o evF{{vLA{i)<yi'ViA{i)) / E) where 
A = —cEeit. Here is the unit vector in the z-th di- 
rection. The conductivity in the i-th direction in the 
lowest energy level approximation valid for t — > oo and 
Vs < Vf, Vd/h is then given by 



a-ii (t) : 



—AevF 



d^kKe[{n"j,{t)\a.\n^^,m+{t)] (9) 



BZ 



with 



f/t't'THt') 



(10) 



=i^[tf dt"- f dt' f dA r"' (t") 

and the transition matrix element 7~"° (t) = 
(u"=(t)|crj|u!=i(t)). By inserting ^ in dll the 
term proportional to t cancels in an improved tight- 
binding approximation since it can be written as 
t/32C?^A:9fc,(u!fi(t)|Ji|u!fi(t)) ^ 0, where J, is the 
tight-binding current operator for A = [8|. Here we 
used the fact that the exact tight-binding wave functions 
are smooth at the Brillouin zone boundary. Summing 



-itAf 



the Fourier series X]aj„ ('^n)^"^"* = 
where Ae = ei — e_i we obtain for large times 



/h^v, (t) 



it 



ujn + Ae/h — iS 



(11) 



In the following we calculate the contribution of every en- 
ergy valley to the momentum integral in ^ separately, 
i.e. o-ji(t) = I]„=o..[v'/7r] '?ii(0(2 - Sn,o)- For large times 
one can restrict the fcj,-integrals of Eq. ([9]) to the neigh- 
bourhood of the valley center ky setting immediately 

ky « ky in Tj"' ■ Then we obtain by partial integration 



d^k 

n— th valley 

ie 1 



{2hvFkxY 
[AeY 



(12) 



+ 00 

dk. 



dk. 



Afc„ 



Vf- 



-J di9ki)Sm\d)[F'''{ki)) - F'^iki))] 
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with 



_(l _ ^-^{2v^k+u^)t^ (J3) 



and fc^ = (tan2(t?) + l)i/2r„{Afc^'^e[cos(i?)] + 
A/Cy '^6[— cos(i9)]} where Q{x) is the Heaviside function. 

In the calculation of an via ([9]), the quanti- 
ties V = dxcxp''^^+~^^^ /d and further Cm = 
Ec^„«2m7r«j./dl^lK)(2 - <5„,o)e'""* are relevant where 

|7^|(a;„) are the Fourier components of |7^|(t). More pre- 
cisely, Cm with TO > are the positive components for 
frequencies 27r(TO - l/2)vp/d < w„ < 27r(m + l/2)v*p/d 
under the restriction that a;„ > for to = 0. A straight- 
forward calculation leads with ([7]) for Vs ^ vp to 



Re[C„i] = ^ cos{C^t)B-X{m,m,m,a 
27rm 



(14) 



-sin(C_t) 



X{m, TO, TO + 1, (t) + Ex, 



Im[C„i] = 2__^ a coa{C'^t)B^X(m,m,m + l,a) 

ae{±} 



+ asm{C-t) 
where 
X{nuj,ni,n2,<y) 



d TT 

sm \A- + ni- 



2TTm 



X{m, TO, TO, (t)-|-Ex 



8(2 



Ad^ Bl - [i{2vtt)B+ + a^Y^ 

-B^+B+a'2v:t)-c+txi2v:t) 



2nv*pt 



1-^ 

V 



no — 
2 



(15) 



functions /i , . . . , /a arc calculated from ([12)) for t ^ oo 

as 



Ii{x,ip) = — -sm^{ip)\cos{(p)\, 

, , 2 2 , „ f aix) 
h{x) —a; -f sgn 1 - Re<^ ^-^ 

TT^ TT^ 4 

0,(2:)'" 



(17) 



a(x) 



-^cr[a(2;)log(l + (Te'"(^))+iLi2(cre*'^(^') 
ae{±} 

z{x) = (^a/ a::^ - 10(2;) - 1^ , 



where a{x) = arctan(l/Va;^ ~ !)■ We obtain from (fTB]) 
that the conductivity ct"^ does not depend on Vg , whereas 
(T^y shows a strong ti<,-dependence. For 



0, d'^y is 



reduced to dyy — (5„.or„e 7r/2ft,. Furthermore we find 
for (T^j. a divergence at SL potentials where V G Ntt for 
general velocities. The same thing holds for dyy but here 
we must demand > 0. The origin of these divergences 
comes from the vanishing of Fq. All this and the result 
of the next section leads us to the remarkable fact that 
an infinite large SL graphene sample is an ideal motion 
detector at SL potentials where the first new Dirac point 
emerges, i.e. at V ~ n: There dyy is vanishing for Vs ~ 
and jumps to infinity for Vs ^ 0. 

From ini), (Ull), the divergence of dyy at V" G Ntt has 
its origin in the approximation that we used an infinite 
ballistic time t ^ ti, in calculating the response. This 
is not really valid for a finite system where tb ~ L/vj 
and L is the length of the sample. By repeating the 
discussion below (jlOp but use the energy eg w sHvpik'^ + 
d^kl/MV^f/^ at y e Ntt leads to d% in ^ with a finite 
cut-off at 1/Fo ~ [VtbY^^ ■ In the following we calculate 
from (fTlTl the conductivities (T° at G Ntt in leading 
order in l/tb for th 00. The results are da ~ d^^ with 



The term Ex in (|T4|) stands for the foregoing expressions 
with interchanged B^ <^ B_, C+ <^ — C_ and switched 
sign of V*. We used further the abbreviation A = A^ + 



A_, B = B+ + B. 



We are now able to calculate the conductivity contri- 
bution (jj^ of the n-th energy valley by using ([9|), (fT2|l . 
dH]), leading to 



^/3 F(l/6)r(l/3)Ai;j.tb\'^' 



h 42/3^ r(2/3) V <^ 
"2 V3 F(l/6)r(l/3)Ai>i.iA'/' 



'^'^y 'h 42/30F 



F(2/3) 



V- 



(18) 

Re[Co+Ci]Re[Coo] 



where Cqo = (vs/d) Jq^^' dtCo. By using (fT4)) we obtain 



;r2 



TT 



-a 



<^yy 



h "°2F„ 

1 TT X - 

" <j&{L,R] 



Rc[Cq+Ci]{ Re [Co 



Re[Ci] 
Im[Ci] 



1 



sin(2C)} - h 



Tn^k^-^d 



(16) 



r„Afc"-'"d 
— y- — ,1?: 



Terms containing C„j with to > 2 are neglected here 
which can be justified numerically. The angle 'd^ is given 
by ian{d'^) = (r„dAfc,;;''")/[7r2 - (F„dAA:^''")2]i/2. The 



32sin(yl|) r 
B+B_Ad^ 



Rc[Coo] = 



cos [{B+-B-}-\-cos (b-^ 



(19) 



In Fig. 2 we plot dyy for Vg/vp = 0.1 (left panel) and 
Vs/vp = 0.01 (right panel), as well as for Vs = (horizon- 
tal curves) at various V^-values. The most interesting V- 
values are where for a certain Vs the signal dyy is largest. 
In particular, the signal to background ratio, i.e., dyy di- 
vided by dyy for Vg = 0, should be large. We obtain from 
the figure and ([T5)) as well as (fT5|) that for a finite system 
and Vs <C vp, the SL potential region where V ^ vp/vs 
and V G Ntt gives the best results. We plot this in Fig. 2 
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FIG. 2: (Color online) We show the conductivity ayy orthogo- 
nal to the SL calculated for velocities Vs/vp ~ 0.1 (left panel) 
and Vs/vf = 0.01 (right panel) as a function of the dimen- 
sionless time. We plot the curves for various SL-potentials 

V by using ^ for V ^ Ntt, (HHJ for F e Ntt and ballis- 
tic times {vptt/dY^^ — 1000. The horizontal curves show 
(Tyy for Vs — 0. Insets in both panels show a zoom in of the 
corresponding curves in the main panels where we add for 

V = 3n/2 an inset showing ayy also for larger times. 



for V = [0.3vf/vs]tt. This is chosen so that the curves do 
not show a higher-order u^-Fourier behavior according to 
([T^ . We note that in principle a graphene velocity de- 
tector based on a SL considered here could also attain a 
large signal to background conductivity for small veloc- 
ity differences by using large SL potentials V ^ vp/ Avs- 
This is due to the phase factors in 

Beside the oscillation frequencies ^ 2TTVp/d and ^ 
2Trv* /d we also find from (fT4|l and Fig. 2 a much smaller 
oscillation frequency ~ C± for the conductivity contri- 
bution of the side- valleys becoming relevant only on very 
large time scales. One can show that due to its non-zero 
velocity, the SL transfer additional energy and momen- 
tum to an electron passing its potential steps such that 
the electron velocity oscillates between ±z;i?(l— and 

±vf[1 — ky{vF ± VsY /^{vf T VsY]- Due to this velocity 
difference the electron picks up an additional oscillating 
phase proportional to t represented by the last term in 
([7]) . This leads to the long wave-conductivity oscillations 
shown in Fig. 2. 

To complete our discussion, we finally calculate the 
quasi-particle velocities in the x and y-direction for elec- 
trons in the w"^-^ state where now £'dc = 0. The knowl- valid for vplkl/ky 
edge of these velocities is useful in quantum pumping 
experiments [HIH] . We obtain from ^ and ^ 



Vf{U^i\<Jxx\U^i) 



VF{u'^l\(Tyy\u'^j_\) 



hdkx ' 



hdky 



5e± 
hdk„ 



(20) 

Re[Co +Ci]. 



trast to Vy. As in the non- moving system [24| there is a 
coUimination of the electron motion in a;-direction, i.e., 
\vy \ ^ \vx\ for potentials were V ~ Ntt and momenta ky 
near the central Dirac point. Here we use that |7^| < 1. 



IV. OUTER- VALLEY TRANSPORT 
CONTRIBUTIONS 

Next, we discuss the conductivity contributions of the 
outer-energy valleys where cfg "C 1. We obtain from ([1]) 
that the new Dirac points are located at kyd/2 = [V^ — 
(ttti)^]^/^ where the linearized energy spectrum around 
these points is given by =shvF[aQk'^+T^{ky — ky)^]^^^. 
The effective y- velocity coefficient is now given by r„ = 
[V^ — {7rn)^]/V^ and do = nn/V. This means that the 
outer-valley regime cfo ^ 1 is fulfilled for those valleys 
where irn/V <C 1. 

We obtain now from Sect. II for the space evolu- 
tion operator 1^ of the non-moving system Xxoi^) ~ 
M sin(a(:^ (a;)2(x — xo)/d)/a^^{x) to leading order in ao- 
The corresponding lowest-band eigenfunctions can be 
interpreted by electrons which are fully backscattercd 
close to the potential steps for \es\d/hvF <C ctQ. This 
is just the opposite situation of the inner-valley trans- 
port contributions discussed in Sect. Ill where we got a 
complete transmission through the potential steps. This 
interpretation is even justified by discussing the scatter- 
ing of electrons on a single potential step in the momen- 
tum regime Sq ^ 1- In this regime u^^ can now be 
written as in ([6]) with the substitution of the spinor part 
{skx/\kx\,l) — > {—isky/\ky\,l). For the moving lattice 
we concentrate us in the following on a particle moving 
in a potential ±V in the region — Wst < a; < d/2 — Vst 
We now determine a complete set of functions fulfilling 
the quasi-relativistic Klein-Gordon equation with a po- 
tential V{x) = ±V obeying the zero-boundary conditions 
v^-Vst,t) = v^{-Vst + d*/2,t) = 0. The distance d* has 
a small modification to the distance d for \es\d/hvF ^ So 
determined by ae^d*/d = nn for the n-th energy val- 
ley, i.e. ao = 7rn. These wavcfunctions consist of a su- 
perposition of two Klein-Gordon wave- function solutions. 
The momenta of both Klein-Gordon wave-functions can 
be formally derived from the zero-boundary conditions. 
More concrete the two corresponding momenta are given 
by a particle initial momentum and its reflected momen- 
tum at the boundary. In the quasi-nonrelativistic limit 
,\ vf wc obtain for these mo- 
menta fc^ ± Vs\ky\/vF and —k^ ± Vs\ky\/vF with j G N 
and kl = 2Trj/d* in the potential ^(a;) = ±V . The 
restriction on the quasi-nonrelativistic limit is justified 
for the outer-valley transport contributions in the case 
Vs ^ Vf- This leads to 



This means that similar to the above conductivity con- 
siderations we obtain no time dependence of Vx , in con- 



yi {x t) = _?_e±*('i"f-\/fex+fc?-^)*/?igTi(i/2)KV"F)|fc«|t 
^ Vd* 

X g±^iv,\ky\/v^)ix+v,t) sin[fcJ (x + Vst)] . (21) 
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By using ([2|) and ([3]), the wavefunction u"" is then given 
by 

d ( ^r, M cosfan) sinfofn) , 9 , 



0.2 0,4 0,6 



8 10 12 14 



2sg[fc, 



sg[kyV{x + Wst)]i 



E 



^sg[V(x+V3t)] j 



(22) 



with 



sg[V(x+v,t)l 



d*/2 



2a; 



(23) 

With this wavefunction in hand we are now prepared 
to calculate the conductivities (t.JJ- for the outer-valleys 
7rn/l/ < 1. By using ^ with (HH) and we obtain 
for the conductivities 







yy 



e TT 1 , , , , , . 
J-,^JUYiv.) (24) 



with Y{vs) is given by 



^(^^)= E -y|c.nc,f ('5.,,+2<5,,>,). (25) 

ij>0 



Here we use \ci\ — \cf\ and Ci,Cj in p5)) and (p3)) are 
calculated with d* — ?> d. 

We obtain from (p4|) that the transport contributions 
of the outer-valleys corresponding to ap 1 show 
no time-fluctuations. This is not based on the quasi- 
nonrelativistic approximation used above. We show in 
Fig. 3 dyy for the outer- valleys and various SL poten- 
tials V and velocity fractions Vg/vp- Most pronounced, 
the curves on the right panel show a conductivity peak 
at valley indices where n sa rip. Here np is given 
by rip [\vs/vf\{V'^ ~ (7rn)2)i/2/7r]. This conductiv- 
ity peak is also observed from ([24|) and ((25|) by tak- 
ing into account that in a rough approximation we have 
|cip ~ (<5i,„+„o + (5i,|„_„o|)/2 leading to 



Y{vs) 



1 



4|7i-noP 4|7i + 7ip|^ 



(26) 



All this means that for V > irvp/vg with ^ 1 we ob- 
tain a large conductivity signal where the conductivity 
modification due to the motion of the SL is of similar 
magnitude as the conductivity value of the non-moving 
SL. Similar is of course true for the detection of small ve- 
locity differences Avg where now we have V > ttvf/^Vs 
in order to obtain a large signal to background value. 
By comparing the conductivity values ayy for the inner- 
valleys ([TB)l . Fig. 2 and the outer- valleys ([M]) . Fig. 3 we 
obtain at least for ^ 1 and V 76 Ntt that the outer- 
valley contributions arc dominant. 



i^lO- \ 1,952 ^ 
tD \ 0,05 0,1 




3ir/2 
.-..in 



-\ ....h 




0,2 0,4 0,6 

vJvf 



FIG. 3: (Color online) Left panel: Outer valley-conductivities 
Gyy (|24p of the n-th electron side- valley as a function oivs/vp 
for certain SL potentials. For V — 57r/2 and V = Stt, which 
both consists of two side- valleys, the upper curve corresponds 
to the valley index n = 1 and the lower curve to n = 2. Inset 
shows a zoom in of o-yy for V 



37r/2. Right panel: ayy for 
V — SOtt as a function of the valley index n for certain SL 
velocities Vs/vf- 



Next, we calculate the effective particle velocities for 
electrons in the outer- valley defined in ([20]) . where now 
again Edc = 0. By using (|2T|) and (|22|1 we obtain 



hdkx 



hdk,. 



(27) 



This shows that there is no Uc-correction term in con- 



trast to the inner- valley case ((20|) for Vy. This is caused 
by the fact that in the outer-valley regime electrons are 
approximately fully reflected, and thus the total proba- 
bility of finding an electron between —Vgi and —Vst + d/2 
is conserved. 

The non-trivial dependence of the conductivities on 
the SL- velocity forced us to treat the conductivity contri- 
butions for the inner and outer-valleys separately. This 
separation is no longer necessary when calculating the 
conductivities for the non- moving SL. For this we use 
the full oscillatory wave- function ([3]) with ([2|) and ([9]), 
PT|) . This leads us to the following Vs = 0-conductivities 



2 

X2"°^ 



1 



^yy 



h2aV 



(28) 



Note here that the magnitudes of ap and F„ correspond 
to the outer-valley values discussed above Eq. pTjl for 
n 7^ and to the n — Q values discussed above Eq. (0). 
Similar expressions as in (|28|l were calculated before 
within the dc vector potential gauge A = 0, leading as in 
pristine graphene to a small overall numerical prefactor 
correction to our result (P5|) [Tsj . The disadvantage of the 
calculation in Ref . [l^ lies in the strong dependence of this 
prefactor on the order of taking the zero-temperature, 
zero- frequency, and zero-damping limit. This does not 
happen in our calculation [sj. 
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V. SUMMARY 

Summarizing, we have considered the dc-transport in 
neutral graphene undulated by a unidirectional moving 
superlattice potential with Vg ^vp, Vd/h. While the re- 
sponse along the direction of the SL wave- vector is van- 
ishing, the dependence is dramatic in the orthogonal di- 
rection. In particular we find for potentials where the 
first new Dirac point emerges, i.e., at = tt, that the 
infinite large graphene sample is a perfect motion detec- 
tor. The orthogonal dc-conductivity is vanishing for zero 
velocity and jumps to infinity at non-zero SL-velocity. 
A large conductivity signal with a high signal to back- 
ground ratio is reached for the finite but large graphene 
system when V € Ntt. The time fluctuating contribu- 
tion to the conductivity is largest when V ~ vp/vs- All 
this was derived from the inner-valley contributions to 
the conductivities. 

Next wc have calculated the conductivity contributions 
of the outer- valleys. The conductivity contributions par- 
allel to the SL-wavcvcctor are vanishing. In the orthog- 
onal direction they arc large, time-independent and ex- 
hibit a peak as a function of the valley index. For V ':$> 1 
and V > -KVp/vs the conductivity modifications due to a 
moving SL are of similar magnitude as the conductivity 
values of the stagnant SL. Note that for F ^ 1 the outer- 
valley conductivity contributions are dominate over the 
inner-valley contributions, at least for V 76 Ntt. Finally, 
we have calculated the conductivities of the non-moving 
SL without the need of a separate calculation for the in- 
ner and outer- valleys. 

Due to its intrinsic low-noise level @ , our results could 
be useful for graphene as a nanophysical motion detector 



device, or even for general sensors based on the surface 
acoustic wave technology f25| . 



Appendix A: Solving the Hamilton-Jacobi equation 
Eq. (HD 

Here wc outline the calculation of ([7]) by solving the 
Hamilton-Jacobi equation (O to first order in fc^. This is 
done with the help of a generalized characteristic method 
plj . The solution is based on the one-particle quasi- 
relativistic orbit x{t) in a moving potential V{x + Vgt). 
With the help of this solution, (x, t) is given by the 
action integral 

Si^{x,t) = f dt'\^^£iL= sg[V{xt{t')+v,t')] 

-V{xt{t')+Vst')^+SlM). (Al) 

Here Xt{t') is the particle trajectory with xt{0) = xq, 
xt{t) = X. The particle- momentum is given by pt'(t) = 
dxS^ {xt{t),t) and the quasi-relativistic velocity by 

xtit') = -VFSg[Vixt{t') + v,t')]^jM^= . (A2) 

Wc note now that it is much easier to determine xt(t') by 
solving the set of equations above for small ky , instead of 
solving the second-order quasi-relativistic Newton equa- 
tion. From this we obtain ([7]). 
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